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We present a calculation of the mass of the rj' meson using quenched and dynamical staggered fermions. We 
also discuss the effects of "cooling" , and suggest its use a quantitative tool. 



i. T] Lore 

One of the more curious inhabitants of the 
low energy spectrum of QCD is the r]' meson. 
Deriving part of its mass from its kinship with 
the Goldstone boson pseudoscalars, the ry' owes 
most of its anomalous weight to its connection to 
the topological susceptibility. Specifically, in an 
SU (3) symmetric world, one would write 

22,2 /I \ 

m^, = TOg + mo (1) 

where mg is the mass common to all of the octet 
meson, and mo is peculiar to the rj'. In the chi- 
ral limit, the octet mass mg vanishes like y/rn^, 
while niQ remains finite. Ignoring mixing and 
identifying m| with the average mass-squared of 
the physical tt's, K's and the rj, one derives the 
"experimental" number of mg = 860 MeV for 
Nf = 3. To complete the calculation of the 
ground state hadron spectrum, one would like to 
compute this number from a first principles lat- 
tice simulation, free of assumptions such as the 
validity of the large Nc expansion. 

In the language of quarks, the ry' is understood 
to gain its mass from "hairpin" diagrams, where 
the quarks annihilate into glue and then recon- 
stitute themselves. In the language of mesons, 
this phenomenon can be described by ascribing 
a strength rriQ to this interaction of flavor singlet 
mesons. Iteration of this process leads to a ge- 
ometric series which shifts the pole in the prop- 
agator from m| to m| -I- mg. In the quenched 
approximation, the series terminates with the sec- 
ond term. The quenched rj' propagator then has 
the unphysical form of a sum of a pole and a dou- 
ble pole, the infrared divergences of which lead to 
extra quenched chiral logs. It is therefore of inter- 
est to see if one can distinguish pole and double 



pole behavior in a lattice simulation. 
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Figure 1. Diagrams for the rj' propagator. 

The quantity of interest is the ratio of the dis- 
connected and connected diagrams. 
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1 — loop 
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If we compute the ratio using A'vai valence quarks 
and iVdyn dynamical quarks, with the same quark 
mass for valence and dynamical, then R(t) should 
asymptote to 



R{t) 



dyn 



exp(— iAm)] 



(3) 



where Am — m^ — mg, and Z' and Z are the 
residues for the creation of the singlet and octet 
particles, respectively. For A^dyn = 0, i.e. in 
the quenched approximation, R{t) should never 
asymptote, but instead rises linearly with slope 
2mQ/mg. 

This approach has been taken with quenched 
configurations and valence Wilson fermions in ref. 
10]. Here we extend the analysis to staggered 
fermions, and use quenched and dynamical con- 
figurations. 
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2. Implementation 

This calculation used three ensembles of lat- 
tices of size 16"^ X 32, one quenched and two in- 
cluding dynamical quarks. The dynamical ensem- 
bles were "borrowed" from the Columbia group, 
and have been used previously in our Bk analysis 
1^ , and the calculation of fs by the MILC collab- 
oration These were generated using molec- 
ular dynamics with Nf = 2 staggered fcrmions 
with quark masses mdyn — -01 and .025. The 
quenched lattices were generated on the Cray- 
T3D at the Ohio Supercomputer Center (OSC), 
using 3-subgroup SU (2) updates, with 4 to 1 mix- 
ture of over-relaxed and heatbath. We saved con- 
figurations with 2000 sweeps, and obtained an av- 
erage update time of 5 microseconds per link on 
16 processors. 

Table 1 
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We computed staggered propagators on the 32 
node T3D at OSC and on the 256 node T3D at 
Los Alamos's Advanced Computer Laboratory. 
We ran our fixed size problem on sets of pro- 
cessors ranging from 16 nodes to 128 nodes, and 
found the per node performance to be essentially 
independent of the partition size over this range. 
By hand coding the inner loops in CAM, the Cray 
T3D assembly language, we obtained a sustained 
performance of 45 Mflops per node, including I/O 
time. 

To create and destroy the rj' we used the stag- 
gered flavor singlet operator 0(75 I)Q- This 
is a distance 4 operator, which we made gauge 
invariant by putting in explicit links, and aver- 
aging over the 24 paths across the edges of the 
hypercube. 

To obtain the one-loop contractions of this op- 
erator (the denominator of fig. |l|), we computed 
two types of propagators. One type used a noisy 
source which was nonzero on two timeslices t — Q 
and t = 1. This source was a random phase 



rjx — e'^(^) for each color and each site, such that 
in the average over noise samples, {rjxvl) ~ 
The second type of propagator used a source ob- 
tained by transporting the noise rjx across the di- 
agonal of the hypercube and putting the phase 
appropriate for the 77' operator. Thus the opera- 
tor which created the rj' was the spatial sum over 
a double timeslice of the hypercube operator, plus 
noise terms which vanish on average. We took 
two noise samples per configuration. 

For the two-loop contractions we used source rj 
with U{1) noise at every site on the lattice. Then 
solving (lp + m)(j) = 77 we estimate the propagator 
as Gxy = {fTKpx'Pl)- We note that this estimator 
is orders of magnitude better than the alterna- 
tive Gxy = {4>xril), but it only available when 
the sites x and y are separated by an even dis- 
tance. For this reason we restricted our attention 
to the pseudoscalar operator, and neglected the 
axial vector (which is distance 3). For each color 
we used 16 noise samples on a doubled lattice, or 
an effective number of 96 noise vectors per mass. 

3. Results 
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Figure 2. Quenched R{t). 

Figure || shows the measured ratio R{t) in the 
quenched ensemble at a valence quark mass rriq = 
.01. In principle, the sickness of the quenched ap- 
proximation should manifest itself in an unending 
linear trend in the data. By contrast, the dotted 
curve is the exponential one would expect if there 
were Nf = 4 active dynamical flavors. Clearly 
the data rule out Nf = 4, and are consistent with 
the quenched form, but from the data alone one 
couldn't rule out a nonpathological behavior with 
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Nf = 2 or fewer flavors. Extracting the slope 
from such curves for several quark masses, we ob- 
tain the values for mg plotted in figure |^ 




Figure 3. Quenched ttiq. 

Extrapolating linearly to ruq — and rescaling 
to the most relevant case of Nf — 3 degenerate 
flavors, we obtain 



]{Nf = 3) = (1050 ± 170 MeV)2( 



2GeV 



)^ (4) 



which is compatible with the "experimental" 
number quoted above. 

In the presence of Nf — 2 dynamical fermions, 
we expect R{t) to asymptote to the constant 2. 
Figure ^ shows the result for rridyn — m^ai ~ 
.01. To fit it to the exponential form may be 
putting more weight on the data than it should 
bear, but its behavior does stand in clear contrast 
to the quenched data. If we ignore our theoretical 
prejudice and simply press ahead to repeat the 
linearized analysis as in the quenched case, we 
find the result 



ml{Nf = 3) = (780±50MeV)2( 



2GeV' 



(5) 



On the other hand, if we do the correct job 
and fit to the exponential form, we can extract 
TO,,/ a at two points, mdyn = mvai = -01 and 
'nidyn = rn^ai — .025, finding .488 ± .030 and 
.676 ± .040 for the bare lattice numbers respec- 
tively. Extrapolating to -niq = and rescaling to 
Nf = 3 we end up with 




Figure 4. Dynamical R{t). 



4. Cooling 

In a closely related study, we have "cooled" the 
dynamical (m^ — .01) to count large instantons. 
Dividing the configurations into bins according to 
their cooled topological charge, we compute the 
rj' mass on two subensembles, one with n = 0,1 
or 2 instantons, and one with n = 3,4 or 5 instan- 
tons. Like the authors of ref. |^ we find that TOq 
mass is larger (by a factor 2) in the subensem- 
ble with larger topological charge, in accord with 
the idea that the rj' gets its mass from instantons. 
This result can be understood by looking at the 
pseudoscalar expectation value in particular con- 
figurations. Figure ^ shows what happens to the 
rj' operator as one cools. The correlation evident 
in the plot persists as a function of cooling time, 
even out to the point where all that remains are 
smooth instantons. We conclude that even in the 
hot configuration, the large features in the rj' are 
due to underlying instantons. 
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ml{Nf = 3) = (730± 250MeV)2( 



2GeV 



)2 Figure 5. Effect of 50 cooling steps. 

This observation invites one to explore a 
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method which could speed up this and other cal- 
culations. Any observable we care to compute 
can be written as 



and we can determine their ratio more precisely 
than we can determine either variable alone. 



(O) = {0),ool X ( 



(7) 



where Ocooi indicates observables computed on 
cooled configurations. In general Ocooi are much 
cheaper to compute, and have smaller statistical 
fluctuations. If the hot and cool observables are 
strongly correlated, it may then be the case that 
for a fixed amount of computing one can obtain 
the RHS of eqn. |^ to better precision than the 
LHS. The method amounts to a type of renor- 
malization group transformation, where Ocooi is 
supposed to retain the long distance physics. 

In a pilot study we tried cooling by only two 
steps. Figure ^ is a scatterplot of TOtt and 
for a variety of quark masses in both the hot 
and cold ensembles. Extrapolating to zero quark 
mass, the hot and cool ensembles give compatible 
rho masses, but evidently the cool data are much 
more precise. Further, if we compare at fixed pion 
mass, the cool data are a factor of 3 to 4 cheaper 
in terms of CG iterations. 
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Figure 6. p and tt masses. 

Picking quark masses which are close to each 
other in figure ^ (m^ = .01 hot, and nig = .04 
cool), we can look for the correlations in more 
complicated observables, such as the disconnected 
{rj' {0)j]' (t)) two-point function. Figure shows a 
scatterplot of the data for t = 10. As the disper- 
sions show, both variables are statistically con- 
sistent with zero, i.e. we are well into the noise. 
Still, the two observables are strongly correlated, 
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Figure 7. Hot and cool observables. 



5. Conclusions and Future Directions 

We have shown that the 77' mass can be com- 
puted using staggered fermions While the data 
are too imprecise for definitive conclusions, it is 
at least plausible that the quenched rj' propaga- 
tor is pathological, while the dynamical one is 
well behaved. To compute the study we want to 
repeat the calculation with smeared sources to re- 
duce the effect of higher excitations. We have also 
suggested a style of calculation which may be of 
some use when computing expensive, propagator- 
intensive observables. Our initial choice of two 
large cooling steps is unlikely to be optimal; to 
fully test the idea one should try better smooth- 
ing algorithms, e.g. using an improved action. 
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